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(54) Verfahren zur quantitativen Analyse von Gasvolumina, tnsbesondere von Abgasen aus 
Verbrennungselnrlchtungen oder -anlagen, sowle Elnrlchtungen zur Durchftihrung der 
Verfahren 

(57) Zur quantitativen Analyse von Gasvolumina, 
insbesondere von Abgasen aus Verbrennungseinrich- 
tungen oder -anlagen mittels Emissions- oder Absorpti- 
onsspektrometrie im ultravioletten, sichtbaren und 
irrfraroten Spektralbereich, werden geometrisch defi- 
nierte und reproduzierbar verstelfbare Beobachtungs- 
ebenen festgelegt, die jeweils senkrecht zur 
L&ngsachse eines Abgasstrahls ausgerichtet sind. 

In einer ersten MeBreihe werden eine Anzahl m 
spektraler Messungen durchgefuhrt, wobei die optische 
Achse eines Spektrometers immer in der jeweiltgen 
Beobachtungsebene liegt jedoch von einer Messung 
zur nachsten jeweils urn eine erste Strecke parallel ver- 
setzt wird. 

In einer zweiten MeBreihe werden n Messungen 
durchgefGhrt, wobei die optische Achse wiederum in 
der Beobachtungsebene liegt, und von Messung zu 
Messung jeweils um eine zweite Strecke parallel ver- 
setzt wird. 

Durch die (m + n) Messungen entstehen zwei 
orthogonale Sichtbundelsfitze. die ein Gitter mit (m x n) 
Schnittvolumina bilden. 

Mit HiHe der (m + n) Messungen ergibt jede Mes- 
sung die spektrale Transmission T(v), bzw. die spek- 
trale Strahldichte l(v), die Qber das gesamte 
Gasvolumen im StrahienbGndeJ des Gesichtsfelds des 
Spektrometers integriert 1st. 
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Beschreiixing 

[0001] Die Erfindung betrifft Verfahren zur quantrtativen Analyse beliebiger Gasvolumina, insbesondere der Abgase 
aus Verbrennungseinrichtungen, wie Triebwerken fur DGsenflugzeuge, deren Brennkammern, andere Gasturbinen. 
5 sowie aus Verbrennungsanlagen, wie Chemieanlagen. Kraftwerksfeuerungen oder Deponien. Ganz allgemein betrifft 
das Verfahren jede Art einer quantitatiyen Analyse von Gasvolumen. Ferner betrifft die Erfindung Einrichtungen zur 
DurchfQhrung der Verfahren sowie Modeller- und Auswerteverfahren. 

[0002] In heutiger Zeit bestimmen u.a. zwei wesentiiche Forderung en die Eigenschaften industrieller Verbrennungs- 
anlagen, zu clenen im vorliegenden Fall auch Flugzeug-Turbinen und ahnliche Verbrennungseinrichtungen gerechnet 

io werden: die Maximierung des energetischen Wirkungsgrads der Verbrennung (Okonomische Forderung), und die Mini- 
mi erung unerwQnschter Bestandleile im Abgas (Okologische Forderung). Okologische Aspekte bestimmen u.a. auch 
das Interesse an der quantrtativen Zusammensetzung der bei naturtichen Brand en entstehenden Abgase. Da Verbren- 
nungsprozesse in aller Regel nicht raumlich homogen ablaufen und Inhomogenrtaten Wirkungsgrad und Abgaszusam- 
mensetzung beeinflussen, ist eine raumlich aufgelOste Analyse der Verbrennung und der Abgase erwQnscht 

is [0003] Analysen finden gegenwartig operationell Qber eine Probennahme statt Dabei kommen sowohl statische 
Samm el verfahren, wie Entnahme von Einzelproben beispielsweise mittels evakuierter SammelgefaBe. als auch dyna- 
mische Sammerverfahren, wie kontinuierliche Probenentnahmen beispielsweise mittels Absaugeinrichtungen, zum 
Einsatz. Eine kontinuierliche. raumlich aufgelOste Abgasanalyse bei den meist groBen Abgasstrdmen der betrachteten 
Systeme ist mit dies en Verfahren nicht Oder nur mit sehr groBem technischen und/oder zeitlichen Aufwand m6glich. 

so [0004] Ferner ist a priori nicht sichergestelft, daB sich die Probe in der Zeit zwischen Entnahme und Analyse, bei- 
spielsweise wegen Temperaturanderungen und/oder chemischen Reaktionen nicht verandert Die Analyse erfordert in 
der Regel urrterschiedliche Verfahren fur unterschiedliche Abgasbestandteile und ist hSufig zertaufwendig. Die Sonden 
zur Probennahme kfinnen unerwOnschte Ruckwirkungen auf den VerbrennungsprozeB haben, bis zur Gefahrdung des 
Betriebs der Anlage, beispielsweise bei Brennkammern oder Triebwerken von Dusenturbinen. 

25 [0005] In neuerer Zeit wind im Rahmen von Fbrschungsprojekten experimentell die Eignung von berOhrungslosen 
infrarot-spektrometrischen Verfahren zur Abgasanalyse untersucht. Mit diesen optischen Verfahren werden zwar meh- 
rere Abgasbestandteile gleichzeitig ertaBt, jedoch integral Qber das Volumen, das vom StrahlenbQndel der Optik des 
Gerats Oberdeckt wind. Daher erlauben sie keine raumliche Auf IGsung, wie sie gerade bei der Urttersuchung von Abga- 
sen aus in der EntwicWung befindlichen Gasturbinen erforderlich ist. sondern nur eine eindimensionale Erfassung der 

50 Abgase. 

[0006] Als nachteilig wird bei gegenwartig eingesetzten Verfahren zur Abgasanalyse angesehen: 

1-Probennehmende Verfahren sind langsam, benfltigen fur jede Abgaskomponente ein eigenes Analysengerat. 
beinhalten die Gefahr der Probenverfaischung, erfbrdern bei hoher raumlicher AuflOsung einen groBen techni- 
35 schen und/oder zeitlichen Aufwand und wirken mCglicherweise auf den VerbrennungsprozeB zurOck. 

2.Ber0hrungslose infrarot-spektrometrische Verfahren ermOglichen nur eine eindimensionale Erfassung der 
Abgase und ergeben die Konzentrationen der Abgaskomponenten im Volumen des Oberdecklen StrahlenbQndels 
nur als Mrttelwert. Gerade bei der Entwicklung von Brennkammern und Turbtnen von Dusentriebwerken ist aber die 
40 Kenrttnis der raumlich aufgelosten Verteilung von Temperatur und Konzent ration der Abgasbestandteile notwendig. 
um Aussagen Qber den Wirkungsgrad und die Emissions* ndizes (Abgasmenge bezogen auf verbrannte Treibstoff- 
menge) machen zu kOnnen. Aus der Kenntnis der raumlichen Verteilung kOnnen auch Eigenschaften der und damit 
Entwicklungskriterien fQr die Brennkammer abgeleitet werden. 

45 3.Die fOr die genaue quantitative Bestimmung von Abgasmengen aus industriellen oder naturlichen Verbrennungs- 
prozessen notwendige raumliche Auf IGsung der Messung ist nicht mfiglich (beruhrungslos) oder aus Grunden des 
erforderlichen Aufwands nicht vetretbar. 

[0007] Wesentiiche Elemerrte der Auswertung sind eine Modellierung der von einem Detektor empfangenen Irrfrarot- 
so Strahlung durch LOsung der Strahlungstransportgleichung und eine Anpassung an die gemessenen Spektren durch 
Ausgleichsverfahren. 

[0008] Bei einer Modellierung des Strahlungstransports (Vorwartsmodell) ergibt sich die von einem Detektor an 
einem Ort s empfangene Strahldichte /(v,s) aus der (integralen Form der) Strahlungstransportgleichung: 



/(v.s) = /(v.s 0 )t(v;s.s 0 ) + J <fc' B(v. 7(s0) d ^* :Sf) (1) 
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[0009] Dabei tsl mit s* der Abstand vom Beobachter auf der Sichtfinie (siehe Frg.1). mit B die Planck-Funktion bei der 
Temperatur T t und mit /fv.sj ein Quetlterm bezeichnet. Die monochromatische Transmission T(v) beschreibt die 
Dampfung der Strahlung zwischen zwei Orten s 0 und s r und ist definiert durch: 



( si 



(2) 



w 
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wobei k(v\s) = k(v ,p(s) t T(s)) der druck- und temperaturabhangige Absorptionskoeffizient und n(s) die molekulare 
TeilchencOchte sind. In inhomogenen Medien sind diese GrGBen ortsabhangig. Die Integration in den Gl/n(l) und (2) 
Qber den Weg wird in den gangigen Strahlungstransportmodellen durch die sogenannte Curtis -Godson-Naherung vor- 
genommen, wobei das inhomogene Medium durch eine Reihe von homogenen Schichten mit konstanter Temperatur 
und MolekOl-Konzentrationen, beispielsweise den Mittelwerten rnnerhalb dieser Schichten, ersetzt wird und Qber diese 
Schichten summiertwird. 

1001 0] Im allgemeinen Fall mit mehreren absoibierenden Molekulen ist k(v\s')n(s') durch eine Summe Ober atle Spe- 
zies Zgk g (>/&)n g (s') zu veraligemeinern. Die von einem Spektrometer empfangene Strahlung ist bekanntlich durch Fal- 
tung des monochromattschen Spektrums (Gl. (1) bzw. (2)) mit der instrumentellen Unienprofiltunktion zu berechnen. 
[0011] Nachteilig bei dieser in den gangigen Auswerteverfahren verwendeten Moddlbeschreibung des gemessenen 
Spektrums ist, daB die endfiche Apertur und der endlichen Blickwinkel (Reld of view FOV) realer Spektrometer nicht 
berticksichtigt wird. Zudem sol It e die "utticherweise" fur Strahlungstransportberechnungen verwendete Curtis-Godson 
Naherung Ner nicht verwendet werden, da die Least-Squares Verfahren, wie unten naher ausgefQhrt, im allgemeinen 
wesentlich mehr Herationsschritte benOtigen als bei LOsung der Integrale mit Standardquadraturverfahren. 
[0012] Die Gl.'n (1) und (2) sind (im allgemeinen nicht lineare) Integralgleichungen der allgemeinen Form 
g(v) = / Y[v,f(s)]ds, die zur Inversion, d.h. zur Bestimmung der unbekannten "Quellfunktion" f, beispielsweise des 
Temperatur- oder Gaskonzentrationsprofils, durch Unearisierung. die auf eine FredholnYsche Integralgleichung erster 
Art, g(v) b iK(v,s)f(s)d$) fQhrt und Diskretisierung in ein lineares Gleichungssystem ubergefflhrt werden. In der Regel 

ist ein Spektrum nur fQr eine endliche Zahl von diskreten Spektralelementen v /( / = 1 m. bekannt In den gangigen 

Auswerteverfahren werden die unbekannten Temperatur- und/oder Gaskonzentrationsprofile nur durch ihre Mitterwerte 
beschrieben, d.h. es werden konstante Profile innerhalb des Abgasstrahls angenommen, beispielsweise ein Ein- 
schichtmodel! des Gasvolumens, Oder ein Dreischichtmodell bei Berucksichtigung von Vorder- und Hintergrund. 
[0013] Falls die Zahl der MeBpunkte m grOBer ist als die Zahl n der Unbekannten, liefert die Diskretisierung der Strah- 
lungstransportgieichung (1) oder G1. (2) ein Qberbestimrrrtes Gleichungssystem. Der OWiche Ansatz ist also eine 
Ldsung durch Minimalisierung, 



mm ||g-y(x)|| 2 = min]T (ft - vWi,*)) 



(3) 
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[0014] Dabei ist g - (g^ t ,g m ) = (g(v,) t . . . ,g(v m )) das gemessene (Transmissions- oder Emissions-) Spektrum 
und der m-Vektor y beschretot das modellierte Spektrum, das im allgemeinen in einem nichtlinearen funktionalen 
45 Zusammenhang mit den unbekannten Mittelwerten fur Temperatur- und Gaskonzentrationen (hier reprasentiert durch 
den n-Vektor x) sowie der Wellenzahl z nach Gi.(1) oder (2))steht 

[001 5] G1 .(3) ist ein nichtiineares Least-Squares Problem fQr den unbekannten n-dimensionalen Vektor x, das durch 
Iteratioa wie GauB-Newton- oder Levenberg-Marquardt- Verfahren, gelOst wird. 

[0016] Ein Nachteil des hier dargestelrten gangigen Verfahrens ist die stark vereinfachte Darstellung des unbekann- 
so ten Temperatur- und/oder Konzentrations- Proffls im Gasvolumen. Zudem muB die Ausdehnung der Schicht als 
bekannt vorausgesetzt werden. 

[0017] Aufgabe der Erfindung ist es daher, Analyseverfahren fur Abgase aus Verbrennungsprozessen zu schaffen, 
welche beruhrungslose Messungen erlauben, schnell arberten, viele Abgaskomponenten und deren Temperatur gleich- 
zertig erfassen und eine Ermittlung ihrer raumlichen Verteilung einschlieBen. Ferner sollen Einrichtungen zur DurchfQh- 
55 rung der Verfahren sowie verbesserte Modellier- und Auswerteverfahren angegeben werden. 

[0018] GemSB der Erfindung ist dies bei Verfahren nach dem Oberbegrrff des Anspruchs 1 , 3 und 4 durch die Verfah- 
rensschritte in dem kennzeichnendem Teil der jeweiligen Anspruche gelGst Vorteilhafte Einrichtungen zur DurchfQh- 
rung des Verfahrens sind Gegenstand der AnsprOche 5 bis 7. 
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[001 9] Femer sind gemaB der Erf indung Modellierverfahren mit den Merkmalen in Anspruch 8 und 9 und ein Auswer- 
teverfahren mit den Merkmalen in Anspruch 10 zur DurchfQhrung des Verfahrens geschaffen. DarQber hinaus and vor- 
teilhafte Werterbildungen der erfindungsgemaBen Verfahren Gegenstand der AnsprOche 1 1 bis 18. 
[0020] Bei den erfindungsgemaBen Verfahren zur quantitativen Abgasanalyse wird ein dreidimensionales, computer- 
gestOtztes Model! der Abgaserscheinung, nflmlich des Abgasstrahis einer Brennkammer, Gasturbine Oder Rakete, 
einer AbgassSule in einem Schornstein, einer Abgasfahne aus einem Vulkan oder Schornstein. usw. verwendet. Ferner 
wind zur DurchfQhrung der Abgasanalyseverfahren gemaB der Erf indung eine spektrometrische MeBeinrichtung einge- 
setzt mtttets welcher die ultraviolette, sichtbare oder infrarote Strahiung am Abgas aus verschiedenen Richtungen 
erfaBt wind. DarQber hinaus wird gemaB der Erf indung ein Auswerteverfahren in Form eines (zwei-, drei- oder) multidi- 
mensionales Inversionsverfahrens eingesetzt. 

[0021 ] Bei gemaB der Erfindung vorgesehenen Abgasanalysatoren wird als MeBgerat ein Spektrometer verwendet. 
welches vorzugsweise im irrfraroten Spektralbereich arbeitet, aber je nach Aufgabenstellung auch for den ultravioletten 
Oder sichtbaren Spektralbereich ausgelegt werden kann. Hierbei wird das Spektrometer in bekannter Weise als Emis- 
sionsspektrometer betrieben. d.h. in volIkDmmen passrver Weise; das bedeutet, es empfangt die von den Abgasen 
emrttierte Strahiung, die in der Regel Infrarotstrahlung sein wird, bzw. die bei hohen Abgastemperaturen beispielsweise 
von Raketentriebwerken auch sichtbare und ultraviolette Strahiung sein kann. 

[0022] Oder das Spektrometer wird in bekannter Weise als Absorptionsspektrometer, d.h. es wind in aktiver Anord- 
nung in Verbindung mit einer Strahlungsquelle betrieben, wobei sich zwischen Strahlungsquelle und Spektrometer das 
MeBvolumen mit dem Abgas bef indet, das von der Strahiung durchlaufen wird, die Strahiung dabei teilweise absorbiert 
und vom Spektrometer empfangen wird. 

[0023] In einer Absorptionsanordnung werden breHbandige Strahlungsquellen im ultravioletten, sichtbaren und vor- 
zugsweise infraroten Spektralbereich, beispielsweise Schwarze Strahler oder auch Globar-Strahler verwendet die auf 
die strahiungsaktiven Spektralbereiche der zu untersuchenden Gase abgestimmt sind. Es werden jedoch auch schmal* 
bandige Strahlquellen, wie beispielsweise Laser und Laserdioden verwendet. Es kOnnen alle bekannten Typen von 
Spektrometern eingesetzt werden. die jedoch f Or die Verwendung in rauher Umgebung geeignet sein mQssen, wie bei- 
spielsweise auf TriebwerksprGfstanderi. 

[0024] Besonders geeignet sind Fourier-Spektrometer nach dem Michelson-Aufbau. vor aiiem wenn sie robust aus- 
gelegt sind. wie beispielsweise MIROR Michelson Interferometer mit rotierendem Retroreflektor, o.a.. Das Michelson- 
Interferometer kombiniert groBe spektrale Bandbreiten mit hohem spektralen NachweisvermGgen (eine Weine rausch- 
fiquivalente Strahldichte NESR). 

[0025] Die erforderliche GrOBe des Spektrometers sol It e im allgemein der naturlichen LinienbreHe der zu betrachten- 
den Gase entsprechen, die bei Umgebungsdruck an der Erdoberfiache ca. 0,1cm" 1 oder besser ist; sie wird im Detail 
jedoch an die spezifischen Gegebenheiten im betrachteten Spektralbereich angepaBt. Beispielsweise kann es erfor- 
derlich sein, eine hOhere spektrale AuflOsung zu realisieren, wenn nSmlich die Linien verschiedener vorhandener Gase 
nahe beieinander liegen und getrennt werden mQssen, urn Wechsetwirkungen zu reduzieren oder zu vermeiden. Auch 
kann eine geringere spektrale Aufl&sung ausreichen, wenn die Linien eines Gases weit auseinander liegen und nicht 
von denen anderer Gase gestOrt sind. 

[0026] Besondere Bedeutung f inden die erfindungsgemaBen Verfahren bei VeitorennungsprozeBen oder thermischen 
Reaktionen, bei denen Abgase freigesetzt werden, sowohl in industriellen Prozessen, als auch in der Natur (Vulkane, 
Waldbrande u.a.). Es kOnnen Abgase beliebiger fester, flussiger oder gasfOrmiger Brennstoffe mit den erfindungs- 
gemaBen Verfahren untersucht werden. Die Abgas-Analyse betrifft die Detektion, Klassifizierung und Quantifizierung 
der im beobachteten Volumen enthaltenen Gase, so wie unverbrannter Brennstoffanteile. wie beispielsweise Kohlen- 
wasserstoffe und teste Partikel wie beispielsweise RuB, sowiedie Bestimmung der Temperatur der Gase und sonstigen 
Bestandteile. 

[0027] Wenn nachfolgend vorwiegend der Begriff Abgase verwendet wird, so sollen damit audh alle anderen Gase 
ausdrucklich eingeschlossen sein, obwohl die Untersuchung von Abgasen als die z.Zt wichtigste Anwendung gesehen 
wind. 

[0028] GemdB der Erfindung ist bei einer Einrichtung zur Durchfuhrung des Verfahrens die optische Achse in Bezug 
auf das Abgas geometrisch definiert und reproduzierbar zu verstellen, so daB damit Beobachtungsebenen durch das 
Abgas zu definieren sind. Diese sind dadurch gegeben, daB mrttels der Einrichtung in einer ersten MeBreihe eine 
Anzahl m spektraler Messungen durchgefuhrt werden, wobei die optische Achse des Spektrometers tmmer in der 
Beobachtungsebene liegt, jedoch von einer Messung zur nachsten urn eine Strecke s x parallel versetzt wird. In ein©- 
zweiten MeBreihe werden n Messungen durchgefuhrt, wobei die optische Achse wiederum in der Beobachtungsebene 
liegt, diesmal aber senkrecht zum Achsenverlauf beim ersten Verfahrensschritt. Hiert>ei wird die optische Achse von 
Messung zu Messung urn eine Strecke s y parallel versetzt. 

[0029] Die Parameter m, n, s x , s y hangen von der Geometrie des Abgasvolumens und der GrOBe des Gesichtsfekfe 
des Spektrometers ab, es kdnnen in speziellen Fallen m = n und s x = s y sein. Die Abstande s x und s y werden so 
gewahlt, daB sie gleich dem Durchmesser des Spektrometer-Gesichtsfeldes am MeBort sind; dann haben die panalle- 
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len SichtbQndel keinen Abstand voneinander, sie kOnnen aber auch Heiner sein, wobei sich dartn die parallelen Sicht- 
bQndel uberiappen, Oder auch grOBer, wobei dann die parallelen SichtbQndel Abstande voneinander haben. Es entsteht 
so ein Satz von (m + n) Messungen. 

[0030] Je nach eingesetzlem Verfahren ergibt jede Messung die spektralen Transmission T(v ) % wobei v[cn>l] die Wel- 
5 lenzahl, d.h. die inverse Weileniange der Strahlung ist, bzw. die spektrale Strahfcfichte l(v), die uber das gesamte Gas- 
volumen im Strahlenbundel des Gesichtsfelds des Spektrometers integriert ist. 

[0031] Durch die Anordnung der optischen Achse in der Beobachtungsebene entstehen bei den Messungen zwei 
orthogonale SichtbQndelsatze, die ein Gitter mit (m * n) Schnittvolumina bilden. Die GroBe dieser Volumina ist durch 
die Geometrie der Anordnung gegeben. wie Raumwinkel des Spektrometers, Entfernung zum Sctmittvolumen. Jedes 

w Schnittvolumen wird dabei aus zwei zueinander orthogonalen Richtungen spektrometrisch erfaBt. 

[0032] Die integrale Messung Qber jedes der m SichtbQndel der ersten MeBreihe enthait n Volumina, die bei der 
orthogonalen zweiten MeBreihe jeweils einem Volumen in jedem deren n Strahlenbundel entspricht Das gilt ebenso 
umgekehrt: Die integrale Messung Ober jedes der n SichtbQndel der zweiten MeBreihe enthait m Volumina, die bei der 
orthogonalen ersten MeBreihe jeweils einem Volumen in jedem deren m StrahlenbQndel entspricht Der Aufbau der 

is MeBanordnung ertaubt es ( analog zur beschriebenen Vorgehensweise weitere Beobachtungsebenen durch das Abgas 
zu legen und in gteicher Weise durch spektrale Messung weitere Datensatze zu gewinnen. Dadurch wird das Abgas 
meBtechnisch in den drei Raumdimensionen erfaBt. 

[0033] AuBer der MeBstrategie mit parallel versetzten optischen Achsen und orthogonalen MeBreihen kann auch eine 
and ere MeBgeometrie gewahtt werden. Beispielsweise kOnnen die MeBreihen mit gefacherten optischen Achsen der 
so SichtbQndel gestartet werden, wobei auch mehr als zwei MeBreihen (Facher) aufgenommen werden kftnnen. Die (nicht 
orthogonalen) Facher fiegen dabei in ein und derselben Beobachtungsebene. Auch hier entstehen wieder Schnittvolu- 
mina, die aus zwei oder mehr Beobachtungsrichtungen erfaBt werden und in bereits beschriebener Weise geometrisch 
bestimmt sind. 

[0034] Die emndungsgemaBen (zwei-, drei- oder) multidimensionalen Auswerteverfahren zur Analyse von Abgasen 
25 verwenden sowohl das eingangs beschriebene Modell der Abgaserscheinung als auch die durch die Messung erhatte- 
nen Datensatze. Es hat zum Ziel. mittels einer Inversionsrechnung die MeBdaten und ihre analytische Beschretoung 
(Ergebnis der ModeOrechnung) in Obereinstimmung zu bringen. wobei sich die gesuchten Parameter (Gaskonzentra- 
tionen, Temperatur, eventuell Druck) ergeben. Dazu bilden die Auswerteverfahren eine raumliche Matrix von Gasvolu- 
mina. deren geometrische Dimensionen (Ausdehnung und Anordnung im Raum) durch die Geometrie der Messungen 
30 und das Spektrometer-Gesichtsfeld bestimmt sind. 

[0035] In der sogenannten Vorwartsrechnung werden mittels des Modells fur jedes Gasvolumen die spektrale Strahl- 
dichte oder der spektrale Absorptionsgrad berechnet, wozu fur die Gaskonzentrationen und die Temperaturen, eventu- 
ell auch den Druck, Schatzwerte angenommen werden, und unter Einbeziehung der Strahlungstransfereigenschaften 
eines jeden Gasvoiumens im SichtbQndel zum Modellwert fur jedes SichtbOndel zusammengefaGt In der Inversions- 
35 rechnung des Auswerteverfahrens werden die modellierten MeBwerte mit den wahren MeBwerten vergiichen und die 
Schatzwerte der zu bestimmenden Parameter (Gaskonzentrationen, Temperaturen, Druck) bis zu einer durch ein 
Abbruchkrrteriumdefinierten optimalen Obereinstimmung iterativ verbessert. 

[0036] Die so gefundenen Parameter stellen anen Mittelwert fur jedes Volumenelement dar, wodurch die 
gewQnschte, raumlich aufgelOste Analyse erreicht ist. Diese wird moglich durch die orthogonale (oder geometrisch 
40 anders strukturierte) MeBstrategie und das (zwei-, drei- oder) multidimensionale Auswerteverfahren. Dabei ist ins- 
besondere vorteilhaft, daB in dem Auswerteverfahren jedes betrachtete MeBvolumen in zweifacher (orthogonaler), 
oder vietfacher (muttidimensionaler) MeBbeobachtung verarbeitet wird. 

[0037] Eine zusatzliche Verknupfung bzw. Parametrisierung der MeBvolumina einer zwei-, drei- oder multidimensio- 
nalen Messung im Auswerteverfahren wird mOglich bei gerichteten und strukturierten Abgasstrfimen. Solche liegen 
45 innerhalb von Schornsteinen sowie bei Dusentriebwerken und Gasturbinen in der Nahe des Dusenaustritts vor. Sowohl 
quer zum Gasstrom als auch in Strdmungsrichtung liegen hier analytisch beschreibbare Abhangigkeiten der Gaskon- 
zentrationen und Temperaturen untereinander vor, die beim Auswerteverfahren genutzt werden kOnnen. Derartige 
Abgasverhaitnisse sind f Or das Analyseverfahren besonders gut geeignet. 

[0038] Anhand der Darstellung einer bevorzugten Ausfuhrungsform werden waiter unten Hinweise fur die Dimensio- 
so nierung eines Spektrometers fur die Abgasanalyse gegeben. 

[0039] Nachfolgend wird die Erfindung anhand von bevorzugten Ausfuhrungsform en unter Bezugnahme auf die anlie- 
genden Zetchnungen im einzelnen eriautert. Es zeigen: 

Ftg.1 schematisch die G^jmetrie der Strahlungstransportrechnung durch ein Medium; 

55 

Fig.2 eine schematische Darsteilung einer erfindungsgemaBen spektrometrisdien MeBanordnung am 

Abgasstrahl eines Flugzeugtriebwerks tm Prufstand und eine mehrdimensionale Strahlungserfas- 
sung; 
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Fig.3 eine schernatische Darsteiiung einer erfindungsgemaBen spektrometrichen MeBanordnung am 

Abgasstrahl eines R ugz eugtr i ebwerks im PrOfstand und eine mehrdimensionale Strahlungserfas- 
sung; 

Rg.4a und 4b einen L&ngsschnitt bzw. Querschnitt durch eine MeBanordnung zur Messung innerhato von Dusen- 



triebwerken im Rug Oder von Gasturbinen im Betrieb; 

Rg.5 schematisch die Geometrie des Strahlurtgstransports bei endlichem Blickwinkel (FOV) des Spektro- 

meters; 

Fig.6 schematisch die Geometrie des Strahlungstransports bei endlichem Blickwinkel und AperturCffnung 

des Spektrometers; 

Rg.7 ein radiales Profil fOr Temperatur- und Teilchendichten der MolekOle; 

Rg 8 eine schernatische Darsteiiung einer Abtastung eines zyTindersymmetrischen Mediums, und 

Rg.9 schematisch eine Abtastung eines allgemein geforrrrten Mediums. 



[0040] Eine MeBanordnung in Rg. 2 zeigt ein Spektrometer 21 in einer bestimmten Entfemung von ein em MeBobjekt, 
in Rg.2 von einem Abgasstrahl 22 eines Rugzeugtriebwerks 23. Der Abgasstrahl 22 ist konzentrisch von einer in Rg.2 
nur schematisch angedeuteten Halterung 24 umgeben, welche in nicht naher dargesteilter Weise so montiert ist, daB 
ste parallel zu sich selbst und konzentrisch zur Langsachse 210 des Abgasstrahl s 22 verschiebbar ist 
[0041] An der Halterung 24 ist eine Anzahl von (m + n) optischen Empfangseinrichtungen 25 montiert, welche nicht 
naher dargestellte, bekannte Linsen- Oder Spiegelsysteme sind. An die (m + n) Empfangseinrichtungen 25 ist je ein 
IJchtlertfaserbGndel 26 gekoppelt, welches an seinem anderen Ende zu je einem Eingang eines optischen Multiplexers 
27 fOhrt Der optische Multiplexer 27 ist mit seinem optischen Ausgang an den optischen Eingang 28 des Spektrome- 
ters 21 gekoppeft 

[0042] Als Lichtieitfaserbundel werden beispielsweise bekannte Quarzfaserbundel im sichtbaren und nahen infraro- 
ten Spektraibereich, Oder Schwermetaflfluorid-Glaser im Bereich 3 fim bis 5 sowie andere fQr cfese und andere 
Spektralbereiche verfOgbare Glaser verwendet Das optische System aus Empfangseinrichtungen 25. Lichtleitfaser- 
bQndeln 26, Multiplexer 27 und Spektrometer 21 ist so dimensioniert daB die Empfangseinrichtungen 25 einen 
Gesichtsfeidwinkel geringer Divergenz (einzelne mrad bis einige Grad) aufweisen. In Abhangigkeit von der geo- 
metrtschen Ausdehnung des MeBobjekts ist der Gesichtsfeldwinkel so dimensioniert, daB Strahlung aus mGglichst par- 
allelen (d.h. mit geringer Divergenz) sich nicht oder nur wenig Gberiappenden StrahlenbQndeln 21 1 errpfangen wird, 
wie beispielhaft in Rg.2 dargestelrt ist. 

[0043] Hierbei sind die (m + n) Empfangseinrichtungen 25 in der Halterung 24 in einer senkrecht zur Langsachse 210 
des Abgasstrahls 22 liegenden Ebene so angeordnet und zueinander justiert, daB jeweils parallel zueinander orien- 
tierte optische Achsen der von den m Empfangseinrichtungen 25 aufgenommenen StrahlenbOndel 21 1 senkrecht zu 
ebenfalis parallel zueinander ausgerichteten Achsen der von den n Empfangseinrichtungen 25 aufgenommenen Strah- 
lenbOndel ausgerichtet sind. So mit spannen die optischen Achsen der (m + n) Empfangseinrichtungen 25 ein aus 
rechtwinWig zueinander veriaufenden StrahlenbQndeln gebildetes Gitter senkrecht zur Achse 211 des Abgasstrahles 
22 auf. In Rg.2 ist zur Vereirtfachung der Darsteiiung m = n = 4 gewShlt 

[0044] Im MeBbetri* (bei laufendem Triebwerk) emHtieren die heiBen Abgase des Triebwerks Infrarotatrahlung, die 
von den (m + n) Empfangseinrichtungen 25 erfaBt, und uber die UchtleitfaserbOndel 26 an den Multiplexer 27 geleitet 
wird. Hierbei erfaBt jede Empfangseinrichtung 25 nur Strahlung aus dem ihr zugeordneten StrahlbOndel 21 1 . Der Mul- 
tiplexer 27 leitet die Strahlung aus alien Faserbundeln sequentiell an das Spektrometer 28 weiter, wo das zugeh6rige 
Spektrum bestimmt wird, d.h. es werden (m + n) Spektren zu jeder MeBposrtion errtlang der Langsachse des Abgas- 
strahls 22 bestimmt. 

[0045] Diese Rohspektren werden in kalibrierte Spektren der Strahldichte bzw. des Absorptionsgrads QbergefQhrt. 
Diese Spektren enthalten spektrale Information aus dem gesamten ihnen jeweils zugeordneten Strahl- oder SichtbGn- 
del 21 1 Dabei stehen die Informationen der einzelnen Bundel untereinander Qber das oben genannte Gitter zueinan- 
der in Beziehung. Jedem der (m + n) StrahlenbOndel ist ein Gasvolumen zugeordnet, dessen Strahlung einmal uber 
dieses Bundel insgesamt integriert, erfaBt und auBerdem in m TeSvolumina errtlang der Bundelachse aufgeteilt wird 
und aus jedem dieser m Teilvolumina enen Strah!ungs(teil)beitrag zu jedem der m orthogonaien StrahlenbQndeln lie- 
fert, und umgekehrt. Diese VerknOpfung der Strahlungsinformation wird in dem anschlieBend beschriebenen Auswer- 
teverfahren so benutzt, wie nachfolgend im einzelnen ausgefOhrt ist. 

[0046] Eine MeBanordnung in Rg. 3 zeigt ein Spektrometer 31 in einer bestimmten Entfernung von einem MeBobjekt, 
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in Fig.3 speziell von einem Abgasstrahl 32 eines Flugzeugtriebwerks 33. Das Spektrometer 31 ist um eine Drehachse 
34 schwenkbar und kann in festen Winkelpoatonen gerasiet werden. Das SichtbQndel 35 wird so in verschiedene 
Richtungen geienkt, in Fig.3 in die Richtungen 351 bis 354. Die Schwenkbewegungen und die Rasterung wird der 
jeweiiigen MeBaufgabe angepaBt, d.h. er kann enger Oder weiter geffichert werden. so daf3 die SichtbQndel 35 sich 
5 Qbertappen kOnnen, luckenlos aufeinander fdgen, Oder LOcken einstellbarer Weite haben; demerrtsprechend ist audi 
ihre Zahl variabel. 

[0047] Eine MeBanordnung in Hg. 4, mit welcher innerhatb von DOsentriebwerken im Rug Oder von Gasturbinen im 
Betrieb gemessen wird, zeigt eine bevorzugte AusfQhrungsform einer MeBvorrichtung gemaB der Erf indung, die fest in 
einem Duserrtriebwerk oder einer Gasturbine instaliiert ist Ein nur schematisch angedeutetes Triebwerk 4 (eine Gas- 

10 turbine) weist ein Gehause (eine auBere Ummantelung) 41 , einen BypaB (Kaltgas) 42, eine HeiBgasduse 43 und einen 
Mischer 44 auf. Die MeBvorrichtung weist Strahlungssender 46, beispielsweise in Form von Schwarzen Strahlern, lich- 
temittiernden Dioden und Einrichtungen 47 zum Strahlungsempfang auf, wie beispieisweise eine Empfangsoptik zur 
Einkoppelung in Licrrtlertfasern, Strah I u ngsd etektor en , u.a. Die Einrichtungen 47 sind radialsymmetrisch so montiert 
und equidistant am inneren Umfang des Gehduses 41 so angeordnet daB stch jeweils ein Sender 46 und ein Empf&n- 

15 ger 47 gegenflberstehen und uber eine Sichtlinie 48 Strahlung vom Sender 46 zum Empfanger 47 geiangt. Die Sichtli- 
nien 48 verlaufen dabei radial zum Abgasstrahl 45. 

[0048] In einer nicht naher dargestelften AusfQhrungsform sind die Sichtlinien der Sender 46 und Empfanger 47 nicht 
radial zum Abgasstrahl orientiert, sondern bilden zueinander parallele MeBstrecken, die senkrecht auf der L&ngsachse 
des Abgasstrahls 45 liegen. Das zentrale Paar aus einem Sender 46 und einer Empfangseinrichtung 47 bildet dabei die 
20 langste MeBstrecke, die gleich dem irmeren Durchmesser des Gehauses 41 ist; tie dazu parallelen MeBstrecken der 
sich anschlieBenden Paare werden vom Zentrum nach auBen stettg kQrzer. 

[0049] Die vorstehend beschriebenen beiden Anordnungen zur Messung in einem Triebwerk (einer Gasturbine) sind 
fur den Betrieb zur Absorptionsrnessung dargestellt. In einer anderen bevorzugten Betriebsweise werden sie auch als 
EmissionsmeB- Anordnungen betrieben. Dazu werden die Sender 46 abgeschattet Oder entfernt Die Empfanger 47 
25 empfangen dann die Strahlung der Gasschichten im Bypass 42 und im HeiBgas 43. 

[0050] In beiden Betriebsweisen werden Messung©! bei ausgeschaftetem Triebwerk (Gasturbine) als Referenzmes- 
sungen und bei eingeschattetem als Betriebsmessungen vorgenommen. Die Verarbeitung von Referenz- und Betriebs- 
messung erfolgt in bekanrrter Weise. 

[0051] ErfindungsgemaB sind Modellierung und Auswertung an die vorstehend beschriebenen MeBanordnungen 

so gemaB der Erfindung angepaBt. 

[0052] Die Modellierung ermOglicht die Berechnung der Strahlungseigenschaften (spektrale Strahklichte, Strahlungs- 
transfer. spektraler Emissions- bzw. Absorptionsgrad) der Abgasbestandteile in AbhSngigkeit von deren rSumlicher 
Dichteverteilung und Temperatur in jeder Volumeneinheit oder Volumenzelle des Abgases, sowie den Strahlungstrans- 
fer durch alle Gasvoiumina bis zu einem Spektrometer. Es modeiliert dabei auch nicht zum Abgas gehOrende Gas- 

35 schichten in der Sichtlinie des Spektrometers (Vorderg-und. Hirrtergrund). Die im Abgas enthaltenen festen und 
gasfOrmigen Komponenten werden zweckmdBig in bekanrrter Weise aus dem verwendeten Brennstoff und dem Ver- 
broinungsprozeB abgeleitet. 

[0053] Fur die Berechnung der Strahlungseigenschaften der Abgaskorrponerrten werden die bekanrtten molekul- 
spektrometrischen Parameter der interessierenden Gase verwendet Unienposrtion, Unienstdrke, Lirienbreite. u.a. 

40 Diese Parameter sind in computerlesbaren Datenbanken wie HIT RAN u.a. gespeichert. Fur den Fall, daB Gase zu 
untersuchen sind, deren molekulspektrometrische Parameter noch nicht bekannt sind, mussen diese in bekanrrter 
Weise in Labormessungen bestimmt werden und kOrtnen dann in dem Verfahren verwendet werden. Das Modell wird 
sowohl zur AbschStzung d^ Strahlungseigenschaften beliebiger Abgasgemische, ats auch im erfindungsgemaBen 
Auswerteverfahren verwendet wie spater gezeigt wird. 

45 [0054] Bei der erfindungsgemaBen Modellierung der beobachteten Strahlung werden die Eirrflusse des Spektrome- 
ters berucksichtigt. Aufgrund des endlichen Offnungswinkels der Speklrometeroptik (FOV. in Rg. 5) sieht dieses nicht 
nur eine einzige Sichtlinie. sondern ein Strahlenbundei endlicher Ausdehnung. Diese Qberiagerung mehrerer Sichtli- 
nien wird durch ein Faftungsirrtegral gemdB 



55 beschrieben, wobei l(v t s. a) die unter dsm Winkel a einfaliende Strahlung ist und die Winkelabhangigkeit durch die 
Funktion FOV(a) beschrieben wird, die ubiicherweise eine Baxfunktion ist Die Integration in Gl.(4) erfolgt durch ein 
geeignetes GauB-Quadratur- Verfahren, wie das GauB-Legendre-Quadratur-Verfahren mit mindestens drei Sichtlinien. 
Wegen der endlichen Empfangsapertur des Spektrometers sind zudem nicht nur ein Strahlenbundei, sondern mehrere 
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parallel versetzte BQndel entsprechend Fig. 6 zu berucksichtigen, was in analoger Weise durch Integration Qber die 
Apertur behandelt wind. (Bei Absorptionsmessungen ist die Modeilierung der Transmission (GL2) in analogs Weise zu 
modifizieren.) 

[0055] Zur Vermeidung der mit der Curtis-Godson-Nfiherung verbundenen Nachteile erfolgt die Integration in den 
5 Gln.(1) und (2) gemaB der Erfhdung durch Standard-Quadraturverfahren, beispielsweise Newton-Cotes Fbrmeln. 
[0056] Bei den erfindungsgemaGen Auswertung werden einerseits mehrere Spektren gletchzeitig ausgewertet und 
artdererseits die Modeilierung der unbekaraiten Profile verbessert. Bei dem erfindungsgemaBen MeBverfahren wend en 
mehrere Spektren gemessen. die im Sinne eines "Global Fit" in einem Least-Squares Verfahren ausgewertet werden. 

Hierzu werden die den Einz el spektren errtsprechenden Me3vektoren g 8 fur s = 1 S zu einem Gesamtvektor 

io G = (g t ,~. ,g 9 ) konkateniert und die Minimalisierungsbedingung wind auf diesen Vektor angewandt: 

min||G-Y(x)|| 2 = min^f;^ - y.^.x)) 2 . (5) 
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(Der Vektor x bezeichnet die zu bestimmenden unbekannten Parameter. Der Einfachheit halber wurden hier Spektren 
gleicher Lange vorausgesetzt.) 

20 [0057] in den derzeit gebrfiuchlichen Auswerteverfahren wird der inhomogene Abgasstrahl durch eine homogene 
Schicht konstanter Bedingungen angenflhert und es werden charakteristische GrOGen dieses Einschichtmodells, z.B. 
die mrttiere Temperatur, CO, C0 2 , oder NO Konzentration bestimmt Bestandteif der Erfindung ist die Bestimmung der 
rflumlichen Verteilung dieser GrOBen. FOr die Diskretisierung der unbekannten Queltfunktion werden erfindungsgem- 
mSB drei Method en verwendet, die in den n&chsten Abschnitten nfiher eriautert werden: 

25 

Parametrtsierung der Queltfunktion durch einen "einfachen" analytischen Ausdruck; 

Diskretisierung durch Errtwicklung nach einem Funktionensystem; 

30 • Diskretisierung durch Quadraturverfahren. 

[0058] Das erfindungsgemSBe Verfahren zur Bestimmung von Konzentrations- und Temperaturprofilen aus den 
gemessenen StraWdichtespektren wird in den Fallen vereinfacht. in denen zusatzliche Informationen Qber die Profile 
vorhanden sind. So treten in der Praxis Falie auf, in denen aufgrund der Kenntnis der chemischen-Zphysikalischen Vor- 
35 gdnge bei der Ausbreitung der Abgase bekannt ist, daB die sich ergebenden Profile durch eine bestimmte Klasse von 
Funktionen beschrieben werden kOnnen. 

[0059] Die Vorgehensweise in solchen Fallen wird nachstehend anhand des Beispiels von Flugzeugstrahftriebwerken 
ertautert. Bei diesen Triebwerken wird das Abgas mit hoher Geschwindigkeit emittiert Eine Durchmischung mit Umge- 
bungsluft f indet daher erst allmahlich statt. Der Ausbrertungsvorgang in den erst en Metern hinter der DOse ist deswe- 
40 gen in erster Linie vom Zustand des Gases beim Austritt abhangig und laBt sich mit Hilfe von Gleichungen beschreiben. 
Es existieren Simidationsprogramme, mit denen Temperatur- und Konzentratonsprof ile basierend auf der Kenntnis der 
Bedingungen in der Austrittsebene vorhergesagt werden kOnnen. 

[0060] Triebwerke kOnnen in Klassen eingeteilt werden, die sich durch die Form der Terrperatur- und Konzerrtrations- 
profile unterscheiden, welche in den ersten Metern hinter der Duse auftreten. Diese Tatsache ist bekannt FOr ziviie 
45 Triebwerke ergeben sich beispielsweise radiale Profile, wie sie in Fig.7 dargestellt sind. Wie der Darstellung des radia- 
len Profits in Fig. 7 zu entnehmen ist sind in der M'rtte (um r=0) die Bedingungen konstant Dieser Bereich verkfeinert 
sich mit zunehmendem Abstand von der Duse bis er ganz verschwindet. Hiefbei fallen weiter nach auBen hin die Gr6- 
Ben entsprechend einer GauB-Funktion ab. 

[0061] Die radialen Profile tassen sich formelmaBig mit Hilfe von Gl.(6) beschreiben: 



{/mu fiir |r| < r c 
, . / 1\t\-tA*\ . . (6) 
U max ~* /u 



fumg) exp (- (|r| r 2 rc)2 ) + U |r| > r c 
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[0062] Dabei ist mit / die Abgastemperatur bzw. die Anzahtdichte der Gase (Anzahl der MolekGle pro Volumeneinheit) 
bezeichnet. Die Indizes 'max' und *umg' Kennzeichnen die Maximalwerte bzw. die Umgebungsbedingungen. r ist die 
Radiustoofdinate, wobei r = 0 im Zentrum der DQse liegt. r c tst der Radius, bis zu dem die Gaseigenschaften konstarrt 
sind. Fur Radien r > r c falK das Prof rl gauBfOrmig mrt dem Halbwertsradius r h auf Umgebungsbedingungen ab. 

5 [0063] Die Parameter r c und r h sind Funktionen des Abstandes von der Duse. Mit zunehmender Entfernung wird r c 
Weiner und strebt gegen null. Gleichzeitig nimrrrtdie Halbwertsbrerte (r^J zu. 1m allgemeinen sind die Geometrieparam- 
ter (r c und r h ) auBerdem von der GrOBe / abhangig. Somit kOnnen sich die konkreten Profile je nach betrachteter Spe- 
zies unterscheiden. Oft jedoch kann ein gemeinsames Profil zumindest naherungsweise angenommen werden. 
[0064] Die gemessenen Strahldichtespektren werden dazu benutzt, die Parameter der funktionellen Zusammen- 

io hange von G1. (6) zu bestimmen. Dies geschieht vorzugsweise simuttan fur alle Messungen mit Hilfe der GauBschen 
Ausgleichsrechnung (least squares fit). 

[0065] Zunachst werden die zu bestimmenden Parameter zu einem Vektor zusammengefaGt. Wind beispielsweise 
vereinfachend angenommen, dafi die Umgebungsbedingungen bekannt sind und alle Profile gletche Geometrieparam- 
ter aufweisen, dann ergibt sich 
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x = I r c» r h» ^max' C 1,max' C 2.max»*"' C A/,max] t 7 ) 



wobei T max fQr die maximal e Temperatur und Ci jfnax fQrdie maximale Anzahldichte des Gases / steht. Der Vektor x wird 
nun mit Hilfe bekannter numerischer Algorithmen (Ausgleichsverfahren) bestimmt. 
so [0066] Im letzten Abschnrtt wurde die funktionale Ortsabhangigkeit der Temperatur- und Gaskonzentrationsprofile 
expfizit vorausgesetzt und die unbekannten Profile wurden durch einfache Funktionen parameterisiert Fur Gas vol u- 
mina wettgehend unbekannter Struktur wird diese Einschrankung fallengefassen und die unbekannten Profile werden 
allgemeiner nach einem geeignet zu wahlenden Satz von Funktionen entwickett. 

[0067] Im eindimensionalen Fall einer Zylindersymmetrie wird das unbekannte Profs in radialabhangige Funktionen 
25 entwickett: 

n 

/(r)-£xy* y (r). (8) 

30 

[0068] Dabei ist f(r) das zu bestimmende (Temperatur- oder Gas-) Profil in Abhangigkeit vom Abstand zur Symme- 
trieachse, <>j ist ein geeignet zu wahlender Satz von Funktionen, z.B. lineare Splines, "hat-functions", Oder die Eigen- 
funktionen des harmonischen Oszillators (^(O-expt-i 2 ). ^(rj-r expf-r 2 )...), und x J( j = 1,„.,n sind die n unbekannten 

35 EntwicWungskoeffizienten, die im Least-Squares Sinne entsprechend Gl.(3) oder (5) bestimmt werden. Im allgemeinen 
Fall mehrerer zu bestimmender Profile wird jedes dieser Profile entwickett und der im Least-Squares-Verfahren zu 
bestimmende Vektor x ist durch die Konkatenierung der EntwicWungskoeffizienten aller Profile gegeben. 
[0069] Dieser Entwicklungs-Ansatz wird auch fur Faile benutzt, in denen nur ein Spektrum entsprechend einer Sicht- 
linie gemessen wurde; die Radialkoordinate ist dementsprechend durch die Koordinate entlang der Sichtlinie zu erset- 

40 zen. 

[0070] Im zwei- (oder drei-) dimensionalen Fall wird in Verallgemeinerung der eindimensionalen EntwicWung Gl. (8) 
wird die unbekannte Profilfunktion als Summe 



f0) = £ x/I> y (f ) (9) 



dargesteltt, wobei die n EntwicWungskoeffizienten x 1 x n wiederum durch Minimalisierung der Residuennorm (Gl. 

so (3) bzw. (5)) bestimmt werden. Die von den beiden Raumkoordinaten des Ortsvektors r = (x,y) abhangigen Funktio- 
nen werden z.B. als Produktfunktionen 0(x,y) = <t>(x) <P(y) gewahtt, also 

55 kj 

[0071 ] Der dreidimensionale Fait ist analog. 
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[0072] Numerische Quadraturverfahren beruhen auf der Ersetzung eines Integrals 



durch eine Summe 



wobei die Gewichte <d, und Knoten s,- vom Integrationsintervall und vom gewahtten numerischen Verfahren abh&ngen. 
Wind eine numerische Integration n-ter Ordnung betrachtet, so werden die Wegintegrale in Ql.(1) bzw. (2) genahert 
is durch 

n 



und mit g t = gfv ,j , A ff « o> y Kr> und x ; « ^s y ; wild eine Matrixgleichung erhalten: 



0/ = 2" A f*y Oder g = Ax. (12) 

25 /-1 

so daB nun die Werte der unbekannten Prof iff unktion x } f(Sj) die mit Hilfe der Least-Squares Methode (3) bzw. (5) zu 
bestimmenden n unbekannten Parameter sind. 
so [0073] Wenn zur Bestimmung von Konzentrations- und Temperaturprofilen aus Emissionsspektren eines zyiinder- 
symmetrischen Mediums angenommen werden kann, daB das zu urttersuchende Medium radialsymmetrisch ist wie 
im Falle des Abgasstrahls, so wird erfindungsgemaB durch radiales Abtasten gemaB Fig. 8 Information uber den Pro- 
firvertauf gewonnen. 

[0074] Zur Auswertung wird eine schmale Scheibe senkrecht zur Symmetrieachse des Mediums (Ausbrertungsrictv 
35 tung des Abgasstrahls) betrachtet und in r Kreisringe eingeteiit Dabei wird angenommen, daB innerhalb eines entspre- 
chend dunnen Ringes die Temperaturen und Konzentrationen konstant and. Aus dem gemessenen Spektrum g 1t 
dessen Sichtlinie durch den fluBersten Ring veriauft, werden dann die entsprechenden Zustandsparameter x 1 des 
Mediums in diesem Ring mit Hilfe eines Least-Squares-Verfahrens bestimmt 



min|l6i-yi(xj)|| 2 , (13) 



wobei y-iCx-)) das entsprechende Model) fur die Beschreibung des Strahlungstransports durch den Ring ist und von den 
Parameter x-j abh&ngt. Mit diesen nun als bekannt anzusehenden Parametem x 1 werden anschlieBend aus dem Spek- 
trum g 2 . dessen Sichtlinie durch den fiuBersten und den nachsten inneren Ring veriauft, ebertfalls mit einem Least - 
Squares-Verfahren die Parameter x 2 dieses Ringes bestimmt Dieses Verfahren wird nun sukzessiv auf jeden weiteren 
so inneren Ring bis zur Symmetrieachse des Mediums angewendet, 



™n||&-y#(xi f x 2l -.. l x i -i f xy)[| 2 fiir i«l f ... f r. (14 > 
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[0075] Auf diese Weise werden alle Konzentrationen und Temperaturen in AbhSngigkeit des Abstandes zur Symme- 
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trieachse des Mediums bestirnmt Diese sequentieile Auswertung IS8t sich einfach implementieren und ermdglicht eine 
rekursive Bestimmung der Schichten. 

[0076] Afternativ werden im Sinne des vorstehend erifiuterten Global Fit Verfahrens die einzelnen Spektren zu einem 
GesamtmeSvektor 

G-ioUl eW 

kombiniert und die Gesarrrtheit aller zu bestimmenden Parameter 

(Vektor x = (xf.xj xj) T ) 

simultan gefittet. 

[0077] Wenn zur Bestimmung von Konzentrations- und Temperaturprofilen aus Emissionsspektren eines beliebig 
geformten Mediums keine Zyiindersymmetrie des zu untersuchenden Mediums angenommen werden kann, wifd das 
Medium aus mindestens zwei Richtungen beobachtet urn r&urrdiche Informationen uber den Vertauf von Konzentrat- 
ions- und Temperaturprofilen zu erhalten. Fig.9 zeigt eine eintache Anordnung zur Untersuchung des Mediums, be der 
zwei senkrecht zueinander angeordnete Beobachtungsrichtungen betracbtet werden. Hieibei ist wiederum ein Schnitt 
durch das Medium gezeigt. 

[0078] Das Raster in Fig.9 beschreibt eine Einteilung des Mediums in Bereiche. Die Numerierung der Bereiche ist 
willkQrlich. Zu jedem Bereich mit der Nummer i gehcrt ein Satz von Parametern x j( die bestirnmt werden solien. Jeder 
dieser Bereiche wird in diesem Beispief aus zwei senkrecht zueinander etehenden Richtungen beobachtet. d.h. zu 
jedem Bereich gehOren zwei Spektren g k und g h in denen Informationen uber die Parameter x; enthalten seien. Werden 
in einer Richtung m und in der anderen Richtung n Spektren gemessen, wobei in Fig.9 m = n = 7 ist. so werden 
S f = mxn Bereiche und insgesamt S g = m + n Spektren erhalten. 
[0079] Durch AneinanderhSngen der Parametervektoren x ( zu einem Gesamtvektor 

X"C*W x^) 7 

und der Spektren g k zu einem GesamtmeBvektor 

G = (glgl g T s f 

gelangt man wieder zu der Formulierung eines Least-Squares-Problems wie in GL(5). 

[0080] Wird die QrtsauflGsung durch dichteres Abtasten und/oder die Zahl der Unbekannten (z.B. durch eine hChere 
Zahl von EntwicWungskoeffizienten in G1.(8) bzw. (9) oder mehr Quadraturpunkte in Q1.(11)) erhfiht, so werden die 
beschriebenen Least-Squares-Probleme schlechter konditioniert. Dies ist einerseits durch die zunehmende Uberlap- 
pung der Sichtstrahlen des Spektrometers bedingt, da Spektren benachbarter Sichtlinien immer ahnlicher werden und 
darret tendenziell immer weniger neue Information enthalten. Andererseits steigt im Fall eines allgemein geformten 
Mediums die Zahl der gesuchten Parameter quadratisch mit der Zahl der gemessenen Spektren an. 
[0081 ] Urn trotz einer schlechten Kondrtiortierung den maximalen Informationsgehalt und damit eine hobere Ortsauf- 
losung aus den gemessenen Spektren zu erhalten, werden sogenannte Regulartsierungverfahren an wend et bei 
denen man zusatzliche Informationen Ober den Verlauf der Profile in das Least-Squares Problem einbringt. Solche 
zusStziichen informationen werden durch die Forderung von Giattebedingungen an die Profile oder durch die Forde- 
rung der Monotonie der Profile in bestimmten Bereichen eingebracht 

[0082] Bei der Regularisierung nach Phillips und Tlkhonov wind eine Glattebedingung in Form einer quadratischen 
Nebenbedingung in das Least-Squares-Problem (G1. (3) bzw. (5)) eingefuhrt 



min (||G - Y(x)|| 2 + A 2 ||L(x - xo)|| 2 ) (15) 



wobei L eine Einheitsmatrix oder eine Matrix ist die der diskreten Darstellung von Ortsabieitungsoperatoren entspricht 
Xo eine a priori Annahme Cber die gesuchten Prof Be ist, und X der sogenannte Regularisierungsparameter ist der das 
Gewicht der quadratischen Nebenbedingung bestirnmt. Somrt wird nicht nur der Abstand zwischen Messung und 
Modell. sondern auch beispielsweise die Norm (der Ableitung) der Profile minimiert. Dies entspricht anschaulich der 
Forderung der Glattheit des Prof Us. Ist ein a priori Verlauf der gesuchten Profile nicht bekannt so wird xq = 0 gesetzt 
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[0083] Da das Model Y(x) im allgemeinen nicht linear fet, wird Gl.(15) durch Linearisierung tm Sinne eines GauB-New- 
ton-Verfahrens iterativ gelOst. FQr die Wahl des Regulartsierungspararneters X existieren eine Reihe von Verfahren wie 
GCV, discrepancy principle, oder L-Kurven- Verfahren. 

[0084] Bei dem erfindungsgemafien Verfahren zur quarrtrtativen Analyse von insbesondere inhomogenen Gasvolu- 
mina werden Spektrometer im infraroten, sichlbaren oder uftravioletten Spektralbereich verwendet, die eine Einrichtung 
haben, mit der die Strahlung des betrachteten Gasvolumens aus verschiedenen, zueinander in def inierter Orientierung 
stehenden Richtungen, (beispielweise auch senkrecht zueinander) empfangen wird. Es wird eine Anzahl Spektren 
gemessenen, die unterschiedlich sind, aber Ober die MeBgeometrie miteinander verknQpft sind. 
[0085] Bei dem erf indungsgemaBen Verfahren wird diese Verknupfung durch Modellierung derselben bei der Auswer- 
tung genutzt. Das Modell errthart auch die analytische Beschreibung einer StrahlungsQbertragung der Szene. Mit HiHe 
der Auswertung werden die gesuchten Gasparameter, wie rftumliche Verteilung der Konzentrationen und Temperatu- 
ren, durch verknppften Vergleich der gemessenen mit den Modellspektren und Minimierung der Quadrate ihrer Diffe- 
renzen ermrttert 



Bezuaszeichenliste 
[0086] 
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Spektrometer 


22 


Abgasstrahl 


23 


Flugzeugtriebwerk 


24 


Hatter ung 


25 


Empfangseinrichtungen 


26 


Lichtiertfaserbundel 


27 


Multiplexer 


28 


optischer Eingang v. 21 


29 


optische Achsen 


210 


Langsachse 


211 


Strahbundel 


31 


Spektrometer 


32 


Abgasstrahl 


33 


Flugzeugtriebwerk 


34 


Drehachse 


35 


UchtbQndel 


351-334 


Richtungen 


4 


Triebwerk (Gasturbine) 


41 


Gehause 


42 


BypaB 


43 


HeiBgasduse 


44 


Mischer 


46 


Strahfungssender 


47 


Einrichtung zum Strahlungsempfang 


48 


Sichtiinie 



Fatenfanspruche 

1- Verfahren zur quarrtrtativen Analyse von Gasvolumina, insbesondere von Abgasen aus Verbrennungseinrichtun- 
gen oder -anlagen mittels Emissions- oder Absorption sspektrometrie im ultra violett en, sichtbaren und infraroten 
Spektralbereich, dadurch gekennzeichnet, daB 

a) geometrisch definierte und reproduzierbar verstelibare Beobachtungsebenen festgefegt werden, die jeweils 
senkrecht zur Langsachse eines Abgasstrahl s ausgerichtet sind; 

b) in einer ersten MeBreihe eine Anzahl m spektraler Messungen durchgefOhrt werden. wobei die optische 
Achse eines Spektrometers immer in der jeweiiigen Beobachtungsebene liegt, jedoch von einer Messung zur 
nachsten jeweils um eine erste Strecke parallel versetzt wird, und 

c) in einer zweiten MeBreihe n Messungen durchgefOhrt werden. wobei die diesmal senkrecht zur Achse im 
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Schritt b) ausgerichtete optische Achse wiederum in der Beobachtungsebene liegt und von Messung zu Mes- 
sung jeweils urn eine zweite Strecke parallel versetzt wird, wobei 

d) durch das Anordnen der optischen Achse in der Beobachtungsebene bei den (m + n) Messungen zwei 
5 orthogonal e Sichtbundelsfttze entstehen, die ein G'rtter mit (m x n) Schnittvolumina biWen, deren G*6Be durch 

die Geometrie (Raumwinkel des Spektrometers, Entfernung zum Schnittvolumen) gegeben ist, und wobei 

e) mit Hilfe dK Satzes von (m + n) Messungen je nach eingesetztem Verfahren jede Messung die spektrale 
Transmission T(v). wobei vfcm 1 ] die Wellenzahl (inverse Welleniange der Strahlung) ist, bzw. die spektrale 

10 Strahldichte ergibt. die Qber das gesamte Gasvolumen im Strahlenbundel des Gesichtsfelds des Spektro- 

meters integriert ist. 

2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB die Beobachtungsebene, in welch er die (m + n) spek- 
tralen Messungen einer Messung durchgefQhrt werden. von Messung zu Messung senkrecht zur Lflngsachse des 

is Abgasstrahls in Richtung weg von dem zu vermessenden Triebwerk versetzt wird. 

3. Verfahren zur quantitativen Analyse von Gasvolumina, insbesondere von Abgasen aus Verbrennungseinrichtun- 
gen Oder -anlagen mittels Emissions- Oder Absorptionsspektrometrie im ultravioletten, sichtbaren und irrfraroten 
Spektralbereich, dadurch gekennzeichnet, da8 

so 

a) geometrisch def inierte und reproduzierbar verstellbare Beobachtungsebenen festgelegt werden, die jeweils 
senkrecht zur Langsachse eines Abgasstrahls ausgerichtet sind; 

b^ in einer ersten MeBreihe eine Anzahl spektraler Messungen mittels einer schwenkbaren Spektrorneteran- 
25 ordnung durchgefQhrt werden. wobei deren optische Achse immer in einer ersten Beobachtungsebene liegt, 

und 

C2) in einer zweiten und folgenden MeBreihen spektrale Messungen in zweiten und weiteren Beobachtungs- 
ebenen durchgefuhrt werden. wobei von MeBreihe zu MeBreihe die jeweilige Beobachtungsebene um eine 
so vorgegebene Strecke parallel versetzt wird. wobei 

e^ rrut Hilfe der MeBreihen je nach eingesetztem Verfahren jede Messung die spektrale Transmission T(v) t 
wobei v[cm' 1 ] die Wellenzahl (inverse Welleniange der Strahlung) ist bzw. die spektrale Strahldichte I(v) 
ergibt die uber das gesamte Gasvolumen im Strahlenbundel des Gesichtsfelds des Spektrometers integriert 
35 ist 

4. Verfahren zur quantitativen Analyse von Gasvolumina, insbesondere von Abgasen aus Verbrennungseinrichtun- 
gen Oder -anlagen mittels Absorptionsspektrometrie im ultravioletten, sichtbaren und infraroten Spektralbereich, 
dadurch gekennzeichnet, daB 

40 

a^eine Beobachtungsebene senkrecht zur optischen Achse eines Abgasstrahles geometrisch definiert und 
reproduzierbar durch eine Anzahl einander diagonal gegenQberliegender Strahlungssender und Strahlungs- 
aufnahmeeinrichtungen festgelegt ist; 

45 ba)in MeBreihen jeweils eine der Anzahl Strahlungssender und -empfangseinrichtungen entsprechende 

Anzahl Messungen in der Beobachtungsebene durchgefuhrt werden, und 

62) mit Hilfe der Messungen je nach eingesetztem Verfahren die Transmission T(v), wobei v[cm" 1 ] die Wellen- 
zahl (inverse Welleniange der Strahlung) ist. bzw. die spektrale Strahldichte IM ergibt, die uber das gesamte 
so Gasvolumen im Strahlenbundel des Gesichtsfelds des Spektrometers integriert ist 

5. Einrichtung zur DurchfOhrung des Verfahrens nach Anspruch 1 , insbesondere zum Messen der Abgase von Flug- 
zeugtriebwerken (Gasturbinen). bestehend aus 

55 einer einen Abgasstrahl (22) umgebenden Hafterung (24), die parallel zu sich selbst und senkrecht zur Lfings- 

achse (210) des Abgasstrahls (22) verschiebbar ist, 

einer Anzahl von (m + n) optischer Empfangseinrichtungen (25), welche an der Halterung (24) in einer senk- 
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recht zur Langsachse (210) des Abgasstrahls (22) liegenden Ebene so angeordnet und zueinander justiert 
sind, da8 jeweils parallel zueinander orientierte optische Achsen der von den m Empfangseinrichtungen (25) 
aufgenommenen StrahlenbQndeln (211) senkrecht zu ebenfalls parallel zueinander ausgerichteten Achsen 
der von den n Empfangseinrichtungen (25) aufgenommenen StrahlenbQndeln ausgertchtet sind. so daB aus 
den rechtwinHig zueinander verlaufenden StrahlenbQndeln ein senkrecht zur Achse (211) des Abgasstrahles 
(22) aufgespanntes Grtter gebiWet ist, 

einer der Anzahl (m + n) optischen Einrichtungen (25) entsprechenden Anzahl von Uchtleitfaserbundeln (26), 
die mrt einem Ende an den Empfangseinrichtungen (25) angekoppett sind; 

einem optischen Multiplexer (27), zu dessen Eingfingen jeweils das andere Ende der IJchtJeitfaserbundel 
gefQhrt ist und 

einem Spektrometer (21), an das der optische Multiplexer (27) angekoppett ist. 

6. Einrichtung zur DurchfOhrung des Verfahrens nach einem der AnsprOche 2 oder 3, insbesondere zum Messen der 
Abgase von Rugzeugtriebwerken (Gasturbinen) gekennzeichnet durch eine spektrometrische MeBanordnung, 
die zur mehrdimensionalen Strahlungserfassung bezQglich eines Rugzeugtriebwerks parallel verschiebbar oder 
schwenkbar ist. 

7. Einrichtung zur DurchfOhrung des Verfahrens nach Anspruch 4, insbesondere zum Messen der Abgase von Rug- 
zeugtriebwerken, dadurch gekennzeichnet daB Strahlungssender (46) und Einrichtungen (47) zum Aufnehmen 
von Strahlung radialsymmetrisch und fiquidistant am inneren Umfang eines Gehfluses (41) so angebracht sind, 
daB jeweils ein Strahlungssender (46) und ein Empffinger (47) einander gegenuberliegend angeordnet sind, so 
daB Ober eine Sichtiinie (48) Strahlung von jedem der Sender (46) zum jeweils gegenQberliegend angeordneten 
Empf&nger (47) gelangt 

8. Modellierverfahren zur DurchfOhrung des Verfahrens nach einem der AnsprOche 1 bis 7, bei welchem mrt einer 
ModelHerung die Strahlungseigenschaften der Abgasbestandteile in Abhangigkert von deren rfiumiicher Dichtever- 
teilung und Temperatur in jeder Volumeneinheit oder Volumenzelle des Abgases scwie der Strahlungstransfer 
durch afle Gasvolumina bis zu einem Spektrometer berechnet werden. 

9. Modellierverfahren zum Auswerten der MeBergebnisse nach einem der AnsprOche 1 bis 7, bei welchem der endli- 
che Offnungswinkel einer Spektrometeroptik (FOV) durch ein StrahienbOndel endlicher Ausdehnung als Uberlage- 
rung mehrerer Sichtlinien in Form eines Fattungsintegrals gemSB Gl.(4) berucksichtigt wind, und 

cfie Integration in GL(4) durch das GauB-Legendre-Quadratur-Verfahren mrt mindestens drei Sichtlinien erfolgt, 
wobei die endliche Emptangsapertur des Spektrometers durch mehrere parallel versetzte StrahienbOndel in Form 
einer analogen Integration Ober die Apertur berOcksichtigt wird. 

10. Auswerteverfahren nach einem der AnsprOche 1 bis 9, dadurch gekennzeichnet. daB mindestens ein Spektrum 
gemessen wird und in einem Least- Squares- Verfahren Temperatur und/oder die rfiumiiche Dichteverteilung der 
Abgasbestandteile bestimmt werden. 

11. Verfahren nach den Anspruch 10, dadurch gekennzeichnet, daB mehrere Spektren gemessen werden und meh- 
rere Spektren gieichzertig im Sinne eines "Global RT in einem Least-Squares Verfahren ausgewertet werden, 
indem die den Einzelspektren entsprechenden MeBvektoren g s fOr s = 1 , . . S zu einem Gesanrtvektor 
G * (g 1 g 8 ) konkateniert werden und die Minimalisierungsbedingung auf diesen Vektor angewendet wird. 

12. Verfahren nach Anspruch 10 oder 11, dadurch gekennzeichnet daB zur Verbesserung einer Modeltierung von 
unbekannten Profilen eine Parametrisierung der Quellfunktion durch einen "einfachen*" analytischen Ausdruck vor- 
genommen wird, und in Fallen, in denen zusfttzliche Irrformationen Ober die Profile vorhanden sind und diese sich 
analytisch beschreiben lassen, mittels der gemessenen Strahldchtespektren die Parameter der funktionellen 
Zusammenhdnge simultan fur alie Messungen mrt Htrfe der GauBschen Ausgleichsrechnung (least squares fit) 
bestimmt werden, indem die zu bestimmenden Parameter zu einem Vektor zusammengefaBt warden, und der Vek- 
tor mil Hitfe numerischer Aigorithmen (Ausgieichsvertahren) bestimmt wird. 

13. Verfahren nach einem der AnsprOche 1 bis 4 sowie 10 und 1 1, dadurch gekennzeichnet daB zur Verbesserung 
einer Modellierung for Gasvolumina wertgehend unbekannter Struktur eine Diskretisierung durch EntwicWung 
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eines geeignet zu wahlenden Satzes von Funktionen durchgefOhrt wird. indem 

im eindimensionalen Fall der Zylindersymmetrie das unbekannte Profil in radialabhangige Funktionen entwickelt 
wird 

S n 

und im zwei- (oder drei-) dimensionalen Fall die unbekannte Prof iifunktion ats Summe 

10 

n 

/(r)-£x/I> y (r) 

15 

dargestellt wird, 

wobei mit f(r) das zu bestimmende (Temperatur- oder Gas-) Profil, <f^ ein geeignet zu wfthlender Satz von Funktio- 
nen rst, (linear e Splines, "hat -functions", oder die Eigenfunktionen des harmonischen Oszillators (<J>o(r) - exp(-r 2 ), 
< J > i(0~ r exp(-r 2 ) ..) ist, und Xj (mit j = 1,..., n) die n unbekannten EntwicWungskoeffizienten sind, die entsprechend 
20 Gl.(3) und (5) bestimmt werden, 

und wobei im zwei- (oder drei-) dimensionalen Fait die von den beiden Raumkoordinaten des Ortsvektors f - (x,y) 
abhSngigen Funktionen als Produktfunktionen <P(x,y) - <P(x)<P(y) gewahlt werden, so daB gilt: 



14. Verfahren nach Anspruch 10 oder 11, dadurch gekennzeichnet, da8 zur Verbesserung einer ModelNerung von 
30 unbekannten Prof Hen eine Diskretisierung durch Quadraturverfahren durchgefOhrt wird. 

15. Verfahren nach einem der Anspruche 1 bis 4 sowie 12 bis 14, dadurch gekennzeichnet, daB fur Gasvolumina 
weitgehend unbekannter Struktur unbekannte Profile nach einem geeignet zu wahlenden Satz von Funktionen ent- 
wickelt werden, im allgemeinen Fall mehrerer zu besttmmender Profile jedes dieser Profile entwickelt wird, und der 

35 im Least-Squares-Verfahren zu bestimmende Vektor x durch die Konkatenierung der EntwicWungskoeffizienten 
aller Profile bestimmt wird. 

16. Verfahren nach einem der Anspruche 1 bis 4 sowie 10 und 1 1 . dadurch gekennzeichnet, daB in Fallen, in denen 
keine Informationen ober Profile vorliegen, das Medium aus mindestens zwei Richtungen in einer Ebene beobach- 

40 tet wird, wobei die Sichtbundel jeweils in ein und derselben Ebene liegen, urn raumliche Informationen uber den 
Vertauf von Konzentrations- und Temperaturprofilen zu erhalten, wodurch die gemessenen Spektren wegen der 
sich kreuzenden SichtbQndel ein Raster und damit eine Einteiiung des Mediums in Bereiche bilden, wobei zu 
jedem Bereich ein Satz von Parametern gehOrt. die zu einem Gesamtvektor aneinander gehangt werden, und ana- 
log die gemessenen Spektren zu einem GesamtmeBvektor aneinander gehangt werden, und das so errtstehende 
Least-Squares-Problem nach dem Giobal-Frt-Verfahren geldst wird, 

Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, daB aus Emissionsspektren eines zytin- 
dersymmetrischen Mediums Konzentrations- und Temperaturprofile bestimmt werden, indem durch radiates Abta- 
sten (Messen) Information Ober den Profilvetiauf gewonnen und sequentiell von auBen nach innen ausgewertet 
wird. 

Verfahren nach einem der Anspruche 1 bis 4 und 10 bis 17, dadurch gekennzeichnet daB in Fallen, in denen die 
gemessenen Spektren wegen sich Gberlappender Sichtbundel zu schlecht konditionierten Least-Squares-Proble- 
men fuhren, der maximale Informationsgehalt und damK eine hdhere Ortsauflfisung aus den gemessenen Spek- 
tren dadurch erhalten wird. daS Reguiarisierungverfahren anwendet werden, bei denen zusatzliche informationen, 
nfimlich die Forderung von Giattebedingungen an die Profile oder durch die Forderung der Monotonie der Profile 
in bestimmten Bereichen, Ober den Veriauf der Profile in das Least-Squares-Problem etngebracht werden. 
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Fig.9 
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